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Abstract 

Classical plasmas in thermodynamic equilibrium should be neither para- nor diamagnetic due to the 
action of the Lorentz force. Magnetic confinement, however, is based on the observed diamagnetism 
of laboratory plasmas. The apparent paradox is investigated on the basis of the resistive magneto- 
hydrodynamic equations. It is found that, at least in simple plasma configurations, these equations 
do not permit a solution, i. e. the paradox cannot be resolved. It seems that the Lorentz force is 
a test-particle approximation which is not suitable to describe the interaction of moving particles in 
agreement with the conservation of energy. 

Resume 

Les plasmas classiques en equilibre thermodynamique ne devraient etre ni para- ni diamagnetiques 
a cause de Taction de la force de Lorentz. Pourtant le confinement magnetique des plasmas en labo- 
ratoire est fonde sur leur diamagnetisme observe. Le paradoxe est explore sur la base des equations 
de la theorie magneto-hydrodynamique resistive. On trouve que ces equations ne permettent pas une 
solution unique, en tout cas pour des configurations simples; c'est a dire qu'il n'est pas possible de 
resoudre le paradoxe. II semble que la force de Lorentz represente seulement une approximation pour 
des particules singulieres. Si l'on veut observer la conservation d'energie, il faut admettre que cette 
force n'est pas appropriee pour decrire exactement l'interaction des particules en mouvement. 
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I Introduction 

Feynman claims in his lectures [1] that both dia- and paramagnetism are exclusively quantum me- 
chanical effects. His argument is of a very general nature: Since the classical Lorentz force is per- 
pendicular to the velocity of a charged particle, the energy mv 2 /2 of the particle does not depend 
on the magnetic field. Two boxes with the same number of particles and the same temperature must 
consequently contain the same energy even if one of the boxes is placed in a magnetic field. If the 
particles would alter the magnetic field, for example diamagnetically, the magnetic field energy would 
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change and the total energy could not be the same. For this reason neither para- nor diamagnetism 
can arise as long as it is only the Lorentz force which acts upon the particles. 

On the other hand, hot fusion plasmas are governed by the Lorentz force alone, since quantum- 
mechanical effects are negligible. Because of the rotational direction of the gyrating particles these 
plasmas behave clearly diamagnetically which is the reason why they can be confined by magnetic 
fields. In tokamaks, for example, one observes magnetic field changes, when the confined plasma is 
heated by external sources. The change of the toroidal magnetic flux is regularly monitored with a so 
called 'diamagnetic loop' [2]. It allows to determine the energy content of the plasma [3] and yields 
information that can be utilized to control the position of the plasma in the vacuum vessel [4] . 

It has been pointed out that an ideally confined plasma is not in any contact with material walls 
and, therefore, not in complete thermodynamic equilibrium. This seems to explain why diamagnetism 
may occur classically in special circumstances. A. Schliiter [5] quoting N. Bohr shows how the 
diamagnetism disappears in a plasma which is placed in a homogeneous field and surrounded by 
reflecting walls (Fig. 1). 




Figure 1 Gyration of particles in a box with reflecting walls (homogeneous field) 



Each gyrating particle constitutes a magnetic moment which by superposition would result in a 
magnetic field opposing the external field (diamagnetism). There is, however, an additional opposite 
magnetic moment created by the particles which are reflected at the walls such that the net effect is 
zero. This example confirms Feynman's conjecture, but it is not a generally valid demonstration. 

A counter-example may be produced by considering an inhomogeneous field: A straight wire 
carries a current and is surrounded by a toroidal vessel containing a plasma (Fig. 2). 
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Figure 2 Gyration of particles in a box with reflecting walls (inhomogeneous field) 



In this case - in addition to the magnetic moment of the gyrating particles - we have also a current 
density in the plasma volume parallel to the external current which is due to the particle drift in the 
inhomogeneous field. This current will alter the externally applied field so that the field energy is 
changed by the presence of the particles. As a consequence the energetic state depends on whether the 
vessel is being placed in the inhomogeneous field or not, in contrast to Feynman's general conclusion. 

In this paper we analyze the apparent paradox by applying the ideal and resistive magncto- 
hydrodynamic (mhd) equations to simple plasma configurations, but we find ourselves unable to 
remove the contradiction. It turns out that there is an intrinsic inconsistency between Lorentz force, 
Maxwell's equations, and energy conservation. This leads us to the conclusion that the Lorentz 
force is a test-particle approximation which ignores the back- reaction on the field-producing magnet. 
In most instances this is justified to a high degree of accuracy, as the test-particle interacts with 
typically 10 23 field-producing particles. When it comes to a plasma, however, the field produced by 
the gyrating particles cannot be neglected any longer and the Lorentz force turns out to be insufficient 
to describe the interaction of the particles in motion. 

II The Resistive and Ideal Magneto-Hydrodynamic Model of a Plasma 

The mhd-equations are derived from the Boltzmann-equation applied to an ionized gas which is 
subject to the action of the Lorentz force. Derivations are found in many text-books. Quoting from 
[6] we have from the momentum balance of electrons and ions in a fully ionized hydrogen plasma: 

— ► — ' civ 

j x B = Vp + m l n— (1) 

E + vxB= — (jxB) --+VJ (2) 
en \ J en 

Here we have put n = ni = n e and omitted terms of the order of the mass ratio m e /nii. Furthermore, 
heat conduction and viscosity are neglected. For sufficiently slow processes Ampere's law holds in 



the form: 



rotB = ^oj (3) 

Faraday's law of induction is: 



dB 
~dt 

We need the equation of continuity for the particle density: 

dn 

and the equation for the internal energy: 



rotE = -- (4) 



div (nv) — —-5- (5) 



/ dp f + 2 -2 

2~dl 2~ P dlvv = S + r )3 (6) 

with p = n (Ti + T e ) and / = 3 for a plasma with three degrees of freedom. S is a power density 
which is added to the plasma, e.g., by electromagnetic radiation. Equation © does not contain the 
magnetic field which is a consequence of the Lorentz force acting presumably on the particles. Apart 
from the Joule heating term the equation is the same as the one for an ideal gas in the absence of a 
magnetic field. 

If the electron temperature is sufficiently high, the terms in @ and |JB} containing the resistivity 
•q, which accounts for the momentum exchange between electrons and ions, may be dropped. The 
resulting system of equations is the model of ideal mhd being valid on a time-scale short compared 
to the electron-ion collision time. 



Ill A Linear Theta-Pinch Heated by Radiation 

We apply the ideal mhd-equations with 77 = to a linear Theta-Pinch in equilibrium, the external 
field of which is produced by a superconducting coil (Fig. 3). 




Figure 3 Field and pressure distribution in a Theta-Pinch 



The straight field lines are parallel to the z- axis. From (|TJ and follows for the internal and 
external field components: 

p(r> + 3£-£ (7) 

We assume that at time t = the plasma is heated by switching on a radiation source so that the 
pressure is increased. Because of (0 the magnetic field must change and the plasma radius a defined 
by p (a) = may be displaced. As long as the radiation source is sufficiently weak, the kinetic energy 
of the plasma motion is negligible compared to the thermal energy and JjJ still holds during the 
expansion of the plasma. The total magnetic flux inside the coil remains unchanged because of Q 
as the electric field vanishes at the surface of the superconductor: 

a b 

^l rdr+ f ^l r dr = (8) 
dt J dt w 

a 

We insert © into (gj: 

^(rvB i ) + -^r = V (9) 

where v denotes the radial component of the plasma velocity, and integrate from to a. Together 
with |JSJ| we obtain an equation for the change of the external field due to the expansion velocity of 
the plasma edge: 

1 dB e 1 da 2 



B e dt b 2 -a 2 dt 
An equation for the divergence of the velocity field follows from ©: 



(10) 



v dp | dp _ 2S f + 2p d(rv) 
dr dt f f r dr 



by elimination of the time derivative of the pressure with 10 and : 

1 d{rv) = 2fi S- fB e (dB e )/dt 
r 8r fB 2 -(f -2) nop 

Integration from to a results in a second equation for the boundary velocity 

da 2 _ ?4 Mo 5- f(dB 2 )/dt 



(12) 



dt J fB 2 -(f-2) f i oP 



dr (13) 



which yields together with IjlOfl an equation for the change of the external field due to the applied 
heating power: 

^l(b 2 -a 2 + J 2 Jl* \ = J ^° Srdr (14) 

dt y J f-(f-2)^p/B 2 ) J f-(f-2)n p/B 2 1 ' 

The task is now to solve (lllf) with the velocity as given by H12fl inside a moving boundary as described 
by IjlOfl and (|14ll . The boundary conditions are: 

[dp/dr] r=0 = , p(o) = 0, «(0)=0. 

As initial condition we may choose an arbitrary pressure profile p (r, 0). If p (a) = is to hold at all 
times, we must require that the heating source S vanishes at the plasma boundary. For the sake of 
simplicity we choose 

S(r,t) = ap(r,t) (15) 

where a is a constant. 



We introduce dimensionless variables: 



r 2 = x a 2 (r) , < x < 1 , t=^-r 
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B 2 (3 F f-2 an 2 / , , x da 2 



(16) 
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The transformation rules are 
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(17) 



Equations and (|12|l transform into: 



du _ (\ 6 dB e \ 1 d(a 2 g e ) 

9a; ' ^2 B e dr J a 2 B e dr ( ' 

The boundary conditions following from (|18(l and l|16(l are: 

dP(0,T) = /9(0,t)[1-(1-*)/?(0,t)] [1 + ^/3(0, r)] f 1 - M ^ 

(20) 



c?r 1 \ B e dr 



0(1, r) = , u(0, r) =it(l, r) =0 
We combine H18(l and l|19(l into a single equation by taking u as the independent variable instead of 

This is possible, since (|19|) is an ordinary differential equation which does not depend on x explicitly. 
Upon substitution of (|19|l equation l|21|l is an inhomogeneous quasilinear equation of first order which 
has the characteristic system: 

ld " w-a-^Hi+Wi 2SdB - 



(3 dr \ B e dr 



1 du _ (1 5 dB e \ 1 d(a 2 B e ) 

P 



(22) 



u dr \2 B e dr ) a 2 B e dr 
The initial condition on u results from (|19|l and (|20l) : 



u(x,0)= [\- 



2 B e (0) 




0)dx- x / (3(x, 0)dx I (23) 



By solving the system of ordinary differential equations (|22f) wc obtain a relationship between u, (3, 
and r when we eliminate the initial profiles u (x, 0) and /3 (x, 0) with l|23() . It may be inserted into 
(|f 9|l in order to express u and /3 as functions of x and r by further integration. 

It turns out, however, that the solution of (12211 does not satisfy the boundary conditions H20J1 in 
general. In order to demonstrate this we choose for simplicity the special case 6 = corresponding to 
f = 2. It will become obvious that the difficulty remains in the more physical case / = 3 . As an initial 
(3- profile we take f3 (x , 0) = 1 — x 2 . The initial u- profile becomes with (I23|) : u (x, 0) = x (l — x 2 ) /6 
. By elimination of x we have: 



u(x,O) = -0(x,O) y/l-0(x,O) (24) 



Integration of (12211 yields: 



1 -(3 1-/3 (x, 0) 



<?(0) 



u = «(aj, 0)^^1-/3(3;, 0) (1 - <) , 5 = a 2 B e 

Elimination of the initial profiles from (24 - 26) leads to the result: 

g(0) PVT=]3e- T 
6 5 (i_ / 3(i_ e -r))2 



(25) 
(26) 

(27) 



In order to determine the time evolution of g, we differentiate (12711 with respect to x and evaluate it 
at P (1 , t) = with (Unj) substituted: 



c9x 



6e T 
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We may also differentiate <|18|) with respect to x and evaluate it at x = 1 imposing (|20|l : 

d 2 /3 
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Integration with respect to time yields: 
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(29) 



(30) 



Elimination of the slope of /3 at the boundary from (|30[1 and l|28|l yields a differential equation for g 

~d0(x, 0) 



ldg 
gdr 



dx 



(31) 



from which g (r) may be determined. 

If we insert, however, l|31|l into l|19|) and integrate from the axis to the boundary we find with 
(HOJ: 



/3 dx 



8/3 (x, 0) 
dx 



(32) 



Obviously, this integral equation is only satisfied at t = 0. At later times (3 evolves according to (|18|) 
from the initial profile to 0(x, 00) — > 1, so that i)M2l) cannot hold at all times. 

In view of this result we come to the conclusion that the set of equations (1-6) has no solution 
in general which would satisfy the boundary conditions. As the problem of heating a Theta-Pinch 
plasma in equilibrium is physically well posed, it must have a solution in reality. Evidently nature 
"uses equations" which are different from those formulated in (1-6). 



IV A Slowly Diffusing Theta-Pinch Equilibrium 

Inclusion of finite resistivity removes the conservation of flux inside the plasma, but it does not remedy 
the situation. Starting from an equilibrium (JJJ with a function p (r) the plasma should slowly diffuse 
and ultimately fill the entire volume inside the coil. Spitzer [6] gives an expression for the diffusion 
velocity: 

VDr l = ( 33 ) 

which is easily derived from JIJ and J2J under the assumption of the magnetic field staying constant 
in time. He rightly remarks that this condition is only satisfied in the test-particle approximation 
Hop/B 2 — > 0, but, nevertheless, he claims one paragraph below that l|33"|) is of general validity re- 
stricted only by the exclusion of inertial terms. 



Let us assume that there is no heating source S, but field energy is dissipated into internal energy 
by the diffusion process. Equation JSJ remains unchanged as long as the plasma diffusion occurs 
inside a superconducting coil. We derive the total power balance by taking the scalar product of J2J 
with j and by eliminating the triple product with Q ■ 

E-J =rjj 2 + v- Vp (34) 

By comparison with © we find: 

Introducing the Poynting vector with and (0J) we have: 

f9p ■ f + 2 d i vi P v) + ^d i v(ExB)+-^^ = (36) 



2 dt 2 v ' u V / 2/i <9i 

Integration over the plasma volume up to the coil radius using Gauss' theorem yields 



fffdp 1 dBf\ , /•/ 1 as L 2 \ , n 

/ 77^7- + ^ ^ rdr+ / -Mrdr = (37 

7 V2 dt 2 fk) dt J J \2u dt J y ' 

a 

as the surface integrals arising from the divergence terms vanish. For simplicity we have omitted the 
term with the kinetic energy minv 2 /2 as the diffusion velocity is very small compared to the thermal 
speed. Furthermore, in accordance with @ we have neglected the electrostatic field energy arising 
through the electric field component E r = (dpi/ dr) / en, which provides the confinement of the ions. 
This contribution is negligibly small in laboratory plasmas compared to thermal and magnetic field 
energy. The heating term 77 j 2 does not appear in l|37f) explicitly. The increase in internal energy 
must come from a decrease of the field energy as the system is energetically closed by the condition 
E = at the superconducting surface. 
We insert J2J into igj: 

— (rvBi + rrij) + -^r = (38) 

This equation must be solved together with the local power balance (JHJ: 

l(dp +v dp\_ f _f + 2 19M 
2 [dt + dr) 2 P r dr 

The solutions for the pressure and the internal magnetic field must satisfy the force balance J2J). We 
write the velocity as the sum of Spitzer's diffusion velocity l|33|l and a term accounting for the effect 
of a finite pressure which is not included in (|33|l : 

rv = \-rv p = us + Up (40) 

Bi 

With the abbreviation x — r 2 equation (|3*5|l reads: 

3 K^) + i^i = o (4i) 



dx v * ' 2 at 



and becomes together with (Q: 



Us + JT-2 Up )d-x + TfT^8-t +p {-dx- + ^x 1 )- (42) 

Substituting the internal field with (J7J) into l|41() yields a second differential equation for the pressure: 

Up * + i| + 2 P ^ = ^ + i L^ (43) 
ox 2 at ox 2n Q ox Afi dt 



It is quite obvious that l)42|l and l)43|l cannot lead to the same solution. The characteristic system of 
the inhomogeneous first order equation (|42ll is: 



and that of (|4"5)l : 



us + f u p /( f + 2) f ~ p(du s /dx + du p /dx) 

up (B| - 2/iop) du p / fe + (dBj)/(4dt) 1 j 

Suppose a solution p(x, t) exists and is known. One may then calculate the velocity u p (x, t) by 
eliminating the time derivative of the pressure from (021 and (gSJ: 

" jf + 2)8 (fipp u s ) /dx + / (dB*)/(4dt) 

JBl-{f-2)^p dX (4bj 

o 

Now the characteristic systems l|44[) and (I45f) may be integrated. The four resulting families of 
characteristics must lie entirely in the surface p(x, t). This is, however, not possible, unless at least 
us (x , t) =0 which is excluded at finite resistivity. 

The so called 'slowly diffusing equilibrium', which is intuitively expected in the case of finite 
resistivity, is not obtainable from the resistive mhd-equations when the pressure dependent term 
omitted by Spitzer in (|33|l is included. Only in the test-particle case /j,qp/B 2 — ► an approximate 
diffusion velocity may be obtained from (|33|) at constant magnetic field. In general, the momentum 
balance (1,2) and the power balance © are inconsistent in conjunction with Maxwell's equations (3, 
4). Again, we must conclude that nature uses a different set of equations or, more properly speaking, 
at least one of the laws of nature as codified in (1 - 6) must be incomplete. It should be noted that 
the discrepancy encountered cannot be resolved by inclusion of heat conduction. Its dependence on 
temperature is different from that of the resistivity so that a cancellation of terms is generally not 
possible. 



V Discussion and Conclusion 

The mathematical model of resistive mhd (1 - 6) is, of course, not an exact description of reality. It 
leaves out not only a number of well known effects such as heat conduction, viscosity, thermoelectricity, 
gravity, but it neglects also the finite mass of electrons, the difference of electron and ion density, the 
effects of quantum mechanics etc. The idealizations involved are common practice in the mathematical 
modelling of physical reality, but this should lead only to minor deviations of the predictions from the 
observations on basic features. In hydrodynamics similar approximations are made, but the results 
as derivable from the model equations are in reasonable agreement with observations. 

In magneto-hydrodynamics the situation is different according to our analysis: The mathematical 
model does not permit a prediction in principle, as the equations are internally inconsistent and do not 
yield unique solutions. This is not acceptable even for an idealized model and one must find out the 
reason. The derivation of the hierarchy of equations follows the same principles as in hydrodynamics 
so that an inconsistency is not to be expected at first sight. If it arises nevertheless, it must have to 
do with an inconsistency in the basic interaction law between individual particles as it is described 
by Lorentz force and Maxwell's equations. 

The validity of Maxwells equations has been sufficiently confirmed and cannot be put into doubt 
in the present context. This is also true for the energy principle. The correctness of the Lorentz force, 
when applied to individual particles, seems also to be sufficiently verified. It is, however, practically 
impossible in these experiments to measure the back-reaction of the orbiting particles on the field 
producing magnet. We can, therefore, not exclude that a term in the elementary force law, which 
either cancels or is negligible in test-particle experiments, has escaped the attention. At least one 
observation raises doubts: 

The angular momentum vector of a negatively charged particle gyrating in a homogeneous field 
is parallel to the field vector. If many particles in a circular conductor rotate in the same direction, 
they form a current which produces a magnetic moment such that the total field at the center of the 



loop is decreased. In this case the conductor is unstable as it has a tendency to turn around an axis 
which lies in its plane when a perturbation occurs. This indicates that the particles in the conductor 
are in a higher energetic state than they would be, if the current of negative particles would flow 
opposite to their natural direction of gyration. It is, of course, well known that a current loop in a 
magnetic field has in addition to its self-energy an extra potential energy which is the negative scalar 
product of its magnetic moment with the field. 

The Lorentz force, however, does not predict a difference in the energetic state of an individual 
particle regardless whether it rotates clockwise or counter-clockwise in a magnetic field. Suppose a 
charged particle is attached at the periphery of a rotatable disk, similar to 'Feynman's paradox' [1]. 
As the Lorentz force points towards or away from the axis of the disk, no extra work is necessary 
to reverse the direction of rotation so that the particle's energy - in contrast to the particles in a 
current loop - is independent of the sense of rotation. In the Appendix we show explicitly how the 
Lorentz force is at variance with the conservation of energy when a single particle interacts with a 
superconducting magnet. 

The comparison of a gyrating particle with a current loop seems to point to an inconsistency 
which is probably at the root of the discrepancy which we have found when the Lorentz force is 
applied to a plasma. Because of the large amount of particles involved, we need to know the correct 
force law describing their interaction, not only a test-particle approximation in an external field. The 
back-reaction on the field producing magnet, which, at sufficiently high pressure, is the plasma itself, 
cannot be neglected any longer. If the Lorentz force would be complemented by a suitable term 
making the energetic state of a particle dependent on whether it is in a magnetic field or not, the 
equation of the internal energy 10 would be altered and, hopefully, the discrepancy could be removed. 

Appendix 

A charged particle is attached at the periphery of a rotatable disk (Fig. 4) which is turned by a 
motor. 




Figure 4 Interaction of a charged particle with a superconducting magnet 
In a concentric superconducting ring flows a current which produces a magnetic field perpendicular 



to the plane of the disk. The particle's equation of motion is: 

m^=F m + q(E + vxB^ (A.l) 

The first term is the motor force acting on the particle via the mechanical fixture, the second term 
is the Lorentz force. The power balance is obtained by taking the scalar product of (A.l) with the 
velocity: 

m dv 2 _ _ ~ , . „. 

~2H = m + q ( ^ 

An electric field at the position of the particle is present when the current in the ring changes in time. 
We derive its tangential component from the vector potential of the superconductor: 



2tt 

/'o dl f b cos tp dtp 

dt ^ dt J ( & 2 + a 2 -2abcostp)- 



where a and b are the radii of the particle orbit and the superconductor, respectively. 
The particle produces also a vector potential: 

4 = M_!_ (A.4) 
p 4tt \x-x'\ V ' 

and an electric field when it is accelerated: 

p _ dl p _ n q 1 dv 

dt - 4tt \x-x'\ dt (A - &j 

The tangential component of this field is: 

Ho q 1 / dv x dv y \ 

E ™ = 17 W^A V^^ + ^ C0SLp ) (A - 6) 

where x' denotes the position of the particle. As it is compelled to move in tangential direction we 
have: 



E = Mog 1 

PV 47T If - X 1 



dv v , dip' 

— cos (ip - ip ) + v v — sin (ip-ip) 



(A.7) 



At the position of the superconducting ring the tangential field component of the particle is: 



E MQg ^cos(y-yO+tt y ^siii(y>-V) 

47r (a 2 + b 2 -2abcos(p-Lp'))^ 



Integration over the angle p from to 27r yields the loop voltage: 

27T , , , 

Un Q f -if cos a + -4r sin a , , 

/ _dt vat -bda, a = ip -ip' (A.9) 

47r J (a 2 + b 2 -2abcosa)* 



In the limit a <C b one obtains: 



A dt b 

Obviously, the particle works on the superconductor which carries a current: 

(jlq q dv v a 



Its energy is increased or decreased depending on the direction of rotation: 



W ps = -^v v I (A.12) 



The superconductor cannot sustain an electric field and must consequently compensate the applied 
voltage to zero by changing its current: 

dl 



dt 



(A.13) 



L is its coefficient of self-induction which may be calculated by inserting into (A.f 3) the self- induced 
voltage 2tt bdA sip /dt: 



o 7 Mo dl 

2-Kb- — 

4tt dt 



2tt 



cos <pbdip 



(r 2 + b 2 — 2r b cos<p) z 



dt 



r^b 



or: 



fi b 2 



cos ip dip 



(r 2 + b 2 — 2r b cos <p) 2 



(A.14) 



(A.15) 



LO - J r ^ b 

L is proportional to the radius b. The factor of proportionality depends on the cross-section of the 
current ring and diverges logarithmically for a thin filament. 

Inserting now (A. 10) into (A.13) we find the connection between acceleration of the particle and 
current change in the superconductor: 



dJ 
~dt 



Hoqa dv^ 



AbL dt 

Substituting this into (A. 3) we obtain the electric field at the position of the particle: 



E, 



sip 



2tt 

UdtJ (a 



cos ipbdip 



^fioqa dv v [i Q a 
{a 2 +b 2 -2 ab cos p)^ ~ AhL dt Ab 



for a <C b 



(A.16) 



(A.17) 



Inserting this into the power balance (A. 2) of the motor yields finally: 



_ -, mdv 2 , 



(fi qa\ 2 1 dv 2 v 



\ 46 



2L dt 



(A.18) 



The work done by the motor is apparently independent of the direction of rotation. To achieve, 
however, the same kinetic energy of the particle, more energy is necessary when there is no supercon- 
ductor present. The additional term arising from the Lorentz force is independent of the magnetic 
field and vanishes proportional to b~ 3 because of (A.15). 

The interaction energy W ps (A. 12) which depends on the direction of rotation and on the magnetic 
field is not accounted for in (A.18). It corresponds quantitatively to the potential energy — fl ■ B of a 
magnetic moment in a magnetic field: The average current of the particle on its orbit multiplied with 
the enclosed area is qav v /2 = and the magnetic field at the position of the particle as obtainable 
from B = rot A is ^ I/2b. One would expect that the motor must supply or gain this extra work, 
but the Lorentz force does not allow for it. 
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Figure captions 

Figure 1 Gyration of particles in a box with reflecting walls (homogeneous field) 
Figure 2 Gyration of particles in a box with reflecting walls (inhomogeneous field) 
Figure 3 Field and pressure distribution in a Theta-Pinch 
Figure 4 Interaction of a charged particle with a superconducting magnet 



